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FOREWORD 


One  mission  of  the  Defense  Mapping  Agency  (DMA»  requires  the  determination 
of  near  real  time  geodetic  quality  absolute  positions  in  areas  far  removed  from  data 
processing,  stations.  It  has  been  demonstrated  that  a  few  hours  of  observations  from 
the  Global  Positioning  System  (GPS)  constellation  can  provide  static  point  position 
solutions  accurate  to  within  3  to  4  m  when  processing  is  performed  with  an  uncor¬ 
rupted  broadcast  ephemeris  and  two  frequency  data.  When  the  GPS  is  operating 
with  AntiSpoofing  (AS)  and  Selective  Availability  (SA),  these  high-quality  solutions 
are  not  possible  unless  the  equipment  can  accommodate  cryptokeys  and  the  neces.sary 
positioning  algorithms  reside  within  the  receiver.  Since  it  is  more  convenient  to  oper¬ 
ate  without  the  overhead  of  cryptokeys  and  in  an  unclassified  environment,  this 
report  investigates  whether  high-quality  absolute  positions  are  possible  using  uncon¬ 
ventional  techniques  with  unclassified  equipment  in  the  presence  of  AS  and  SA. 

Support  was  provided  by  DMA  under  the  direction  of  Mr.  B.  Roth  and  Mr.  S. 
Malys. 

This  report  has  been  reviewed  by  Dr.  J.  Blanton,  Head,  Space  and  Geodesy 
Branch;  Mr.  T.  Sims,  Head,  Space  Sciences  Branch;  and  Mr.  J.  Sloop,  Head,  Space  and 
Surface  Systems  Division. 
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ABSTRACT 


This  report  addresses  the  problem  of  determining  absolute  positions  ( longitude 
and  latitude)  from  Global  Positioning  System  (GPS)  data  in  the  presence  of  Selective 
Availability  (SA).  The  SA-epsilon  effect  can  be  avoided  if  one  can  wait  until  corrected 
or  post-fitted  satellite  ephemerides  and  clock  estimates  become  available.  However, 
the  SA-dither  effect  cannot  easily  be  eliminated  from  the  ranging  data  in  the  field. 
This  report  considers  both  problems  and  suggests  methods  that  have  the  potential  to 
allow  absolute  position  solutions  to  be  obtained  even  with  dither  corrupting  the 
satellite  data.  The  requirements  and  practicality  of  each  are  reviewed. 
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MOTIVATION 

Geodetic  surveying  of  Fixed  points  with  the  Global  Positioning  System  (GPS)  uses  the 
satellite  constellation  to  relate  the  coordinates  of  a  point  being  surveyed  to  the  World  Geodetic 
System  (WGS84)  coordinate  system.  The  first  half  of  the  link  between  the  point  and  the 
coordinate  system  depends  upon  the  observations  of  pseudoranges  and/or  phase  to  locate  the 
point  with  respect  to  the  constellation  of  satellites.  The  second  half  of  the  link  locates  the 
satellites  in  space  and  time  with  respect  to  the  monitor  stations  whose  positions  are  precisely 
known  in  the  WGS84  coordinate  system.  With  the  Department  of  Defense  Selective  Availability 
(SA)  in  operation,  both  of  these  links  are  purposely  corrupted.  The  SA-dither  introduces  errors 
that  appear  as  random  range  errors  on  the  observations  and  weakens  the  link  between  the  site  to 
be  surveyed  and  the  satellites.  The  SA-epsilon  introduces  errors  into  the  satellite  ephemerides  and 
weakens  the  link  between  the  satellites  and  the  coordinate  system. 

Because  the  removal  of  the  SA  effects  in  the  field  requires  additional  cost  and  operational 
difficulty,  it  is  assumed  that  the  user  would  prefer  to  operate  without  that  overhead.  GPS 
receivers  that  do  not  have  the  capability  to  remove  the  SA  effects  are  less  expensive  to  purchase 
and  can  be  operated  on  the  Lj  frequency  in  an  unclassified  mode.  Under  normal  operations  with 
these  receivers,  the  SA-epsilon  effect  can  be  avoided  if  one  can  wait  to  do  the  data  processing 
until  post-fitted  satellite  ephemerides  become  available  from  a  third  party.  However,  the 
SA-dither  effect  cannot  be  removed  from  the  ranging  data  except  in  a  secure  environment.  This 
report  addresses  the  dither  problem  and  investigates  methods  that  may  allow  absolute  position 
solutions  to  be  obtained  even  with  dither  corrupting  the  satellite  data. 

The  AntiSpoofing  (AS)  encryption  will  not  be  addressed.  AS  prevents  access  to  the 
precise  code  (P-code)  on  both  Lj  and  L2  frequencies.  The  Coarse  Acquisition  (C/A)  code  is 
available  at  all  times,  but  only  on  the  L,  frequency,  and  is  subject  to  the  SA  degradation 
described  above. 


BENEFIT  EXPECTED  FROM  SUCH  A  TECHNIQUE 

When  SA  is  on  and  AS  is  either  on  or  off,  unclassified  civilian  receivers  can  still  be  used 
to  collect  GPS  ranging  data  to  do  absolute  positioning  if 

1.  Single  frequency  C/A  code  is  used. 

2.  SA-epsilon  is  avoided  by  using  a  post  fit  precise  ephemerides. 


1 


NSWCDD/rR-93/309 


3.  The  measurement  system  is  designed  to  be  insensitive  to  SA-dither. 

The  quality  of  the  absolute  position  solutions  so  obtained  was  examined  based  upon  an 
experiment  performed  at  the  Naval  Surface  Warfare  Center,  Dahlgren  Division  (NSWCDD).  It 
cannot  be  expected  and  it  has  not  been  demonstrated  that  solutions  can  be  obtained  that  are 
comparable  to  the  level  of  accuracy  attainable  when  SA  is  not  in  use. 


TECHNIQUES  REVIEWED 

It  has  been  demonstrated  that  collecting  GPS  data  with  two  or  more  antennas  in  the 
relative  positioning  mode  results  in  solutions  for  short  baselines  that  are  accurate  and  repeatable 
to  within  about  1  mm.  The  differencing  of  the  phase  observations  between  the  antennas  makes 
this  method  insensitive  to  errors  due  to  satellite  position,  propagation  effects,  satellite  clock 
variations,  and  SA-dither.  W'hereas  the  insensitivity  is  a  benefit  for  precise  relative  positioning, 
the  method  .makes  it  difficult  to  say  anything  about  the  absolute  location  of  the  baseline.  Relative 
positioning  gives  a  baseline  vector  that  is  precise  in  length  and  direction,  but  is  free  to  translate 
through  large  distances  in  space. 

Processing  large  quantities  of  GPS  observations  obtained  over  short  baselines  may  help 
resolve  the  location  of  the  observer  if  the  dominant  source  of  error  is  random.  However,  this 
method  is  inherently  weak.  What  is  helpful  for  absolute  positioning  is  some  additional 
observation  that  is  correlated,  in  an  absolute  sense,  with  a  point  on  the  earth.  This  auxiliary 
observation,  along  with  the  vector  baseline,  would  be  used  to  pinpoint  the  location  from  which 
the  observations  were  made.  The  acceptable  accuracy  for  this  new  scheme  must  be  better  than 
that  available  from  GPS  single-frequency  Standard  Positioning  Service  (SPS)  when  SA  is  active. 
It  should  be  able  to  routinely  achieve  20  m  Root  Sumed  Square  (RSS)  with  a  goal  of  5  m. 

Geodetic  quality  absolute  positioning  requires  precise  pseudorange  or  range  difference 
data  and  precise  satellite  ephemerides  and  clock  estimates.  The  absolute  position  of  the  unknown 
point  is  found  by  collecting  data  from  several  satellites  over  a  period  of  time.  This  provides  the 
link  from  the  user  to  the  satellites.  The  satellite  positions  are  known  from  the  satellite 
ephemerides  that  are  uploaded  by  the  GPS  Control  Segment  and  transmitted  to  the  user  in  real 
time.  A  globally  distributed  array  of  monitor  sites,  whose  positions  are  known,  receive  data  from 
the  satellites  on  a  continuous  basis.  These  data  allow  the  past  satellite  positions  to  be  estimated 
and  the  future  positions  to  be  predicted.  These  ephemerides  provide  the  link  between  the  earth 
and  the  satellites.  In  effect,  the  user  is  positioning  an  unknown  site  with  respect  to  the  global 
array  of  monitor  sites  with  the  GPS  satellites  acting  as  an  intermediary. 
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When  relative  positioning  techniques  are  used,  the  absolute  link  to  the  monitor  sites  is 
weak  and  an  alternate  link  to  the  earth  is  required.  Several  auxiliary  observations  that  have  been 
considered  include  the  following: 

1.  Physically  position  the  baseline  so  that  it  is  normal  to  the  geoid.  The  normal  to  the  geoid 
could  be  transformed  to  the  equivalent  normal  to  the  ellipsoid  through  the  application  of 
local  deflections  of  the  vertical.  If  this  could  be  accomplished  precisely,  the  observed 
vector  baseline  would  take  on  a  direciicn  that  is  unique  to  some  particular  point  on  earth. 
The  geodetic  longitude  and  latitude  could  then  be  computed  through  the  equation  for  the 
ellipsoid  normal. 

2.  Pointing  the  GPS  baseline  at  some  object  whose  absolute  position  is  known  also  gives  a 
means  for  finding  the  location  of  the  observer,  if  the  range  can  also  be  determined.  The 
range  can  be  found  by  electronic  distance  measurement  devices.  This  scheme  can  be  very 
useful  in  a  local  area  where  some  previously  surveyed  landmark  is  prominent. 

3.  Employing  astronomical  objects  as  angle  references  could  be  another  scheme.  The  GPS 
antennas  mounted  on  a  telescope  or  other  pointing  device  could  be  aimed  directly  at  some 
astronomical  object  whose  coordinates  are  known.  Observations  of  several  baseline  vectors 
directed  at  different  objects  should  be  enough  to  uniquely  determine  the  observer’s  location. 

4.  Precise  optical  astronomical  observations  are  always  dependent  upon  weather  conditions; 
therefore,  a  similar  technique  utilizing  all-weather  radio  observations  would  be  more 
desirable.  Radio  interferometry  may  be  used  to  replace  the  optical  measurements  described 
above  with  the  positions  of  the  GPS  satellites. 

5.  Collect  data  from  one  or  more  baselines  separated  by  moderate  distances  and  perform 
conventional  double  difference  processing  to  determine  the  baseline  vectors,  hold  the 
vectors  fixed,  and  reprocess  the  data  with  the  absolute  position  of  a  reference  antenna  as 
the  unknown.  No  auxiliary  measurements  are  necessary  for  this  technique. 

6.  GPS  and  leveling  data  are  collected  along  two  baselines  radiating  outward  from  a  fixed  site 
and  preferably  separated  by  about  90  deg.  Tangents  to  the  WGS84  ellipsoid  model  are 
formed  by  subtracting  the  leveling  difference  in  height  along  the  baseline  from  the  GPS- 
determined  height  and  correcting  for  both  the  Earth’s  curvature  and  deflections  of  the 
vertical.  The  latitude  and  longitude  of  the  fixed  point  is  found  from  the  ellipsoid  normal 
vector.  The  normal  is  computed  from  the  cross  product  of  the  two  tangent  vectors  obtained 
from  the  baseline  vectors. 

A  brief  description  of  each  technique  will  be  presented  following  the  preliminary  analysis. 
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PRELIMINARY  ANALYSIS 

The  formulation  presented  below  outlines  the  procedure  for  processing  the  data  received 
by  the  GPS  receiver(s)  from  the  two  antennas.  The  observation  equation  is  written  in  the  form 
suggested  by  Braasch:* 

p=p+c(r  -r3  +  d  *  d  +J,+ 

r  s  '■  rcvr  jv'  tropo  lono  rtl  H) 

URE  +  d  „  +  +  noise  +  SA 

mp  nW  RkUS 


where  the  symbols  represent  the  following: 


Pr 

P, 


^  tropo 
drtl 

URE 


mp 

dm 

dfngas 

noise 

SA 


observed  pseudorange  as  observed  by  the  receiver  (m) 
true  geometric  range  between  satellite  and  antenna  (m) 
receiver  clock  offset  from  GPS  time  (sec) 
satellite  clock  offset  from  GPS  time  (sec) 
propagation  delay  caused  by  the  troposphere  (m) 
propagation  delay  caused  by  the  ionosphere  (m) 
satellite  relativity  correction  (m) 

User  Range  Error  (satellite  ephemeris  and  clock  errors  projected  onto  the  range 

direction)  (m) 

multipath  error  (m) 

receiver  hardware  delay  (m) 

receiver  measurement  bias  errors  (m) 

receiver  measurement  noise  (m) 

Selective  Availability  epsilon  and  dither  (m) 


For  the  purpt)ses  of  the  report,  the  following  terms  will  be  ignored:  URE,  d^^^ 

and  noise.  The  SA  term  will  be  expanded  into  its  comptment  parts:  SA^  =  SA-epsiton  and  SAf,  = 
SA-dither.  The  remaining  terms  will  be  retained  but  subscripts  will  be  added  to  indicate  time  i, 
satellite  j,  and  antenna  k. 

~  ^  ^ 1  *  (2) 

dionojk^h^  ^  dmeasjk  *  d„pji^+  ^^(0 


In  Equation  (2),  the  term  is  the  magnitude  of  the  difference  between  the  vector  to 

satellite  j,  represented  by  and  the  vector  to  the  antenna  receiving  the  signal  rjt)- 


Psjk  =  I 


(3) 
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Next,  collect  data  from  the  same  satellites  with  the  two  antennas.  Difference  the  data 
between  the  two  antennas  to  remove  the  satellite  dependent  variables,  including  both  components 
of  SA.  The  differenced  observation  equation  is 


AW  =  -  >‘.12(0  * 


mttujl 


d  d.  -  d-  .. 

mttuj2  tonojl  k>noj2  mpjl 


(4) 


Continuing  with  the  second  difference  between  satellites  m  and  n  results  in  Equation  (5). 
Now  there  is  no  contribution  from  the  local  clock. 


*miasml 


*mtasm2 


measnl 


^iono  ml  ^iono  m2  ^ 

^mp  ml  ^mp  m2 


Id, 

id, 

f  d,„pnl 


*mtasn2 


1  + 


'iononl  d,0„gn2  ■ 

dmpn2  ^ 


(5) 


When  two  receivers  are  used  and  synchronized  by  the  same  frequency  standard,  double 
differencing  to  remove  receiver-dependent  terms  may  not  be  necessary.  In  this  case,  there  is  just 
one  local  time  offset  variable  instead  of  two  and  and  that  term  in  Equation  (4) 

disappears.  What  is  left  is  integer  count  offsets  for  each  of  the  satellites  and  differences  between 
the  measurements  noises,  refraction  effects,  and  multipath.  The  ranges  P.1  and  functions 

of  the  X,  y,  z  components  of  the  satellite  positions  and  the  coordinates  of  the  antennas. 


POLE  NORMAL  TO  THE  GEOID 

One  solution  to  the  problem  of  minimizing  the  effects  of  SA-dither  on  absolute  position 
solutions  requires  that  a  short  baseline,  whose  end  points  are  defined  by  the  phase  centers  of  two 
GPS  antennas,  be  oriented  so  that  it  is  normal  to  the  geoid.  Since  the  phase  centers  of  the 
antennas  are  not  precisely  known  with  respect  to  their  physical  structure,  there  will  be  an 
unavoidable  uncertainty  in  identifying  when  the  baseline  is  in  fact  normal  to  the  geoid.  It  must 
be  normal  with  a  high  degree  of  accuracy  because  it  defines  a  unique  vector  at  that  particular 
location.  This  vector  is  unique  because  there  is  only  one  vector  normal  to  the  geoid  at  any  given 
point.  A  diagram  of  the  situation  is  shown  in  Figure  1.  Tliree  surfaces  are  illustrated:  the 
equipotential  surface  called  the  geoid,  the  mathematical  surface  called  the  ellipsoid,  and  the 
physical  surface  labeled  the  terrain.  The  two  vectors  at  the  point  Pq  represent  the  ellipsoid  normal 
and  the  geoid  normal.  The  angular  difference  between  the  two  normals  is  called  the  deflection 
of  the  vertical  e.  It  can  be  seen  that  e  changes  from  px^int  to  point  due  the  continuously  changing 
directions  between  the  two  normal  vectors. 
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FIGURE  1.  NORMALS  TO  THE  ELLIPSOID  AND  THE  GEOID  AT  A 

POINT  P„ 


Observed 
Geoid  Normal 


True  Geoid 
Normal 


Terrain . 


Ellipsoid' 


True  Ellipsoid 
/  Normal 

When  the  Observed  Geoid  Normal 
is  Corrected  for  the  Deflection 
of  the  Vertical,  It  Becomes  the 
Observed  Ellipsoid  Norrr 


FIGURE  2.  ERROR  0  IN  THE  DETERMINATION  OF  THE  GEOID 

NORMAL 
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In  order  to  determine  the  geodetic  absolute  position,  the  normal  to  the  ellipsoid  at  must 
be  known.  Assuming  that  the  baseline  formed  by  the  GPS  antennas  can  be  made  to  coincide  with 
the  local  normal  to  the  geoid,  GPS  will  provide  the  geodetic  com{X)nents  for  that  vector;  x^,  y^, 
and  Zg  as  in  Equation  (6).  Additional  information  about  the  deflection  of  the  vertical  is  then 
required  to  transform  the  true  geoid  normal  into  a  true  ellip.soid  normal  as  shown  by  the 
application  of  the  angle  e  in  Figure  2.  From  the  approximate  absolute  position  solution  t^btained 
using  the  conventional  GPS  SPS,  the  vertical  deflection  angles  '%  and  q  (where  accurately  known) 
can  be  found  from  a  table.  The  angles  ^  and  q  are  the  east  and  north  components  of  e.  After  this 
operation,  the  ellipsoid  normal  unit  vector  u  can  be  used  with  Equation  (7)  to  solve  for  the 
geodetic  longitude  /.  and  latitude  (j>. 


r,  =  ♦  yp  -  (6) 

u  =  cos^  cosyl  X  +  cos4>  sin/l  y  +  sin^  £  (V) 

If  the  geoid  normal  is  not  fully  determined  and  as  a  result  there  is  an  error  0  in  its 
direction  as  illustrated  in  Figure  2,  application  of  the  corrections  for  the  deflections  of  the  vertical 
will  propagate  the  error  into  the  observed  ellipsoid  normal;  then  the  result  is  an  error  0  in  the 
ellipsoid  normal  vector.  Under  these  conditions,  when  Equation  (2)  is  u.sed  to  determine  the 
longitude  and  latitude  of  the  point  P„,  the  error  may  be  con.>iderable  because  the  position  for 
which  the  observed  ellipsoid  normal  is  the  true  normal  may  be  quite  distant.  This  is  illustrated 
on  an  ellipsoid  with  an  exaggerated  eccentricity  in  Figure  3.  The  place  where  the  observed 
ellipsoid  normal  is  actually  normal  to  the  ellipsoid  is  the  point  P,.  In  addition,  the  residual  errors 
in  the  deflection  of  the  vertical  at  P„  will  degrade  the  accuracy  of  this  process. 

In  order  to  put  the  problem  on  a  firmer  foundation,  a  numerical  example  follows.  Suppose 
that  the  error  in  the  determination  of  the  baseline  vector  is  1  mm  in  the  direction  perpendicular 
to  the  true  normal.  Assuming  a  baseline  length  of  1  m,  the  angular  error  incurred  is  1  mrad 
(angle  0  in  Figure  2).  If  the  deflection  of  the  vertical  is  perfectly  compensated,  then  the  error  in 
the  ellipsoid  normal  is  also  1  mrad.  For  numerical  simplicity,  the  point  Py  will  be  placed  at 
longitude,  latitude,  and  ellipsoidal  height  equal  to  zero.  This  is  on  the  equator  at  the  intersection 
of  the  plane  formed  by  the  x  and  z  axes  of  the  WGS84  coordinate  system.  In  this  case,  the  true 
unit  vector  from  Equation  (7)  is  The  observed  unit  vector,  which  includes  the  1-mrad  error, 
is  r/„  with  components  jc„„=0.9999995,  y^^z^O.OOl,  and  z„„=0. 

u=xx  +  yy-^zz  (^) 

O  UO  ^140  ^  ^ 


Solving  by  equating  components  in  Equations  (7)  and  (8),  gives  the  latitude  as  zero 
because  there  is  no  £  comp.inent.  Solving  for  longitude  from  either  of  the  other  two  components 
gives  an  observed  longitude  of  0.05729578  deg.  This  is  equivalent  to  a  1-mrad  longitude  error. 
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which  is  the  same  magnitude  as  the  error  in  the  ellipsoid  normal.  A  similar  result  in  latitude  is 
obtained  if  all  the  error  is  put  in  No  matter  how  the  error  is  distributed,  the  angular  position 
error  on  the  ellipsoid  is  equivalent  to  the  magnitude  of  the  angular  error  in  the  baseline  vector. 
In  this  example,  if  the  magnitude  of  the  vector  from  the  origin  to  the  pxiint  Pq  is  given  by 
Equation  (9),  then  Vq  is  6378135  m,  and  the  position  error  on  the  surface  of  the  ellip.soid  is 
approximately  6378.135  m. 


(9) 


This  technique  is  very  sensitive  to  errors  incurred  when  the  GPS  antennas  are  set  up  to 
form  a  baseline  perpendicular  to  the  geoid.  A  small  error  in  this  placement  translates  directly  into 
a  position  error  equal  to  that  angular  error  multiplied  by  r„.  The  lact  that  the  ellipsoid  normal  and 
the  observed  ellipsoid  normal  differ  by  a  distance  approximately  equal  to  the  angular  error  0 
multiplied  by  the  height  of  the  terrain  above  the  ellipsoid  (h)  is  not  relevant  when  Equation  (2) 
is  being  used  to  determine  the  geodetic  position.  Even  in  the  case  where  the  terrain  and  the  geoid 
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coincide  with  the  ellipsoid  (h=0),  an  angular  error  6  in  the  observed  normal  vector  would  still 
cause  an  error  in  the  determination  of  the  geodetic  position. 

In  this  application,  with  GPS  baselines  on  the  order  of  1  m,  the  errors  in  the  vector 
direction  will  be  on  the  order  of  1  mrad.  Given  that  the  other  errors  will  be  added  to  this,  this 
scheme  does  not  appear  to  be  practical.  An  alternate  approach  to  the  data  processing  problem  is 
presented  next. 

Rapp^  gives  the  transformation  that  connects  north,  east,  and  vertical  at  a  given  longitude 
and  latitude  with  Earth  Centered  Earth  Fixed  x,  y,  z;  this  is  shown  as  Equation  (10)  below.  The 
and  are  observations  of  the  vector  baseline  components  obtained  from  GPS  and  corrected 
for  the  local  deflections  of  the  vertical.  If  we  admit  that  no  matter  how  careful  we  are  positioning 
the  baseline  normal  to  the  geoid,  there  will  always  be  small  errors  in  the  east  and  north  directions 
whose  components  are  and  These  might  be  accounted  for  if  they  are  treated  as  bias 
parameters  that  must  be  included  in  the  position  solution.  From  each  baseline  solution,  there  are 
three  observations  y*,  and  zj  to  be  used  on  the  left-hand  side  of  Equation  (10).  There  are  five 
unknowns  (X,  <j),  Cj,,  and  on  the  right-hand  side.  Therefore,  at  a  given  site  there  must  be 
k  sets  of  observations  from  different  satellite  geometries  that  give  k  independent  baseline 
determinations  for  the  components  x*,  and  z^.  For  each  baseline  solution,  there  are  3k 
observations  and  5  unknowns.  When  hi2  the  number  of  observations  equals  or  exceeds  the 
number  of  unknowns  and  a  solution  for  longitude,  latitude,  and  the  biases  should  be  possible. 

Xj  =  -(sin^os/1) /ij  -(sin>l)e|,  +(cos(pcosA) 

y^  =  -(sin0sin>l)nj  +(cosyi)e^  +(cos^inA)Mj  (10) 

=  (cos^)nj,  +(sin0)«^ 


An  initial  guess  for  X  and  (j»  can  be  obtained  from  the  SPS  ptisition  solution.  With  the 
baseline  approximately  normal  to  the  ellipsoid,  the  initial  guess  for  and  is  zero  and  for 
is  the  measured  baseline  length.  A  least  squares  solution  including  two  or  more  observations  of 
the  baseline  vector  may  be  able  to  improve  on  the  initial  estimate  of  K  and  (Ji.  Ultimately,  the 
accuracy  of  the  method  will  depend  upon  the  knowledge  of  the  local  vertical  deflections. 


POINTING  AND  RANGING  TO  A  KNOWN  LANDMARK 

If  a  GPS  survey  is  desired  in  a  localized  area  that  has  one  or  more  precisely  known 
absolute  references,  it  would  be  possible  to  use  a  single  receiver  of  the  direction-finding  type 
along  with  an  electronic  distance  measuring  device  to  locate  points  in  line-of-sight  to  the 
references.  A  more  accurate  way  to  perform  the  same  task  (with  two  receivers)  would  be  to 
locate  one  GPS  receiver  over  the  reference  pt^int  and  another  successively  over  the  unknown 
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points.  If  this  were  not  possible  for  some  reason  or  if  only  a  single  direction-finding  receiver 
were  available,  the  survey  could  be  performed,  but  with  lower  accuracy  than  with  two 
independent  receivers.  This  method  is  the  inverse  of  that  described  by  O’Leary,  Evans  and 
Smith. ^ 


Figure  4  illustrates  the  geometry  that  the  method  requires.  The  reference  point  is  located 
by  the  known  WGS84  vector,  r^,  and  is  indicated  by  the  solid  dot.  The  point  to  be  determined 
is  indicated  by  the  +.  The  unit  vector  from  the  unknown  point  to  the  reference  point  is 
determined  by  aiming  a  theodolite-like  device  at  the  reference  point.  The  theodolite  would 
include  a  short  GPS  baseline  positioned  over  the  unknown  point.  The  range  R  to  the  reference 
point  is  obtained  from  a  laser  rangefinder.  Errors  in  the  measurements  are  represented  by  AR  and 
A0  in  Figure  4.  O’Leary  concluded  that  the  angular  error  introduced  1  m  of  error  per  kilometer 
of  range  (2-m  GPS  baseline).  The  range  error  was  found  to  be  about  5  m  for  noncoop)erative 
targets  up  to  a  limiting  range  of  a  few  kilometers.  A  noncooperative  target  implies  that  no  laser 
retroreflector  was  used. 


POINTING  TO  A  STELLAR  OBJECT 

This  technique  depends  entirely  upon  angular  measurements;  therefore,  several  sightings 
to  stellar  objects  are  required  for  good  results.  The  method  is  similar  to  that  exploited  by  the 
sextant,  but  instead  of  measuring  the  elevation  of  an  object  above  the  local  horizon,’’  the  three- 
dimensional  unit  vector  to  that  object  is  determined.  The  theodolite-like  device  equipped  with  a 
short  GPS  baseline  is  aligned  with  the  stellar  object.  Real  time  GPS  attitude  algorithms  can  be 
used  to  obtain  the  instantaneous  unit  vector. 

With  several  unit  vectors  to  known  stellar  objects  found,  a  least  squares  fit  to  the  data  can 
be  used  to  determine  absolute  position.  Assuming  precise  alignment  of  the  optical  axis  with  the 
GPS  baseline,  the  accuracy  of  the  result  would  depend  upon  transverse  baseline  error.  If  that 
error  were  about  1  mm,  then  the  angular  error  for  a  1-m  baseline  would  be  about  1  mrad.  The 
results  would  be  similar  to  those  described  in  the  next  section  with  the  exception  that  attitude¬ 
determining  GPS  receivers  could  be  used  unmodified  and  with  current  algorithms  to  determine 
the  unit  vectors.  Additional  work  would  be  necessary  to  integrate  these  algorithms  with  the  stellar 
data  and  with  the  final  absolute  positioning  routine. 


INTERFEROMETRIC  GPS  OBSERVATIONS 

Standard  GPS  relative  positioning  techniques  are  sometimes  referred  to  as  interferometry. 
Classic  interferometry  implies  that  the  array  of  antennas  form  an  array  reception  pattern  by 
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Reference  Point  — 

Ranae  Error 

AR 

•••  / 

.•••■/  \ 

-  An^le  Error 

.y 

y 

Point  to  be 

Determined 

A0 

HGURE  4.  GPS  AND  RANGE-FINDING  TECHNIQUE 


adding  or  multiplying  the  antenna  signals  together  before  processing.  In  standard  GPS  relative 
positioning,  the  observations  from  each  antenna  are  treated  independently  and  incoherently.  The 
Doppler  phase  counts  are  referenced  to  the  receiver’s  local  frequency  standard  and  the  numeric 
phase  differences  are  recorded  at  regular  intervals.  In  interferometry,  the  received  radio  frequency 
phases  are  added  directly  so  that  they  interfere  and  form  the  array  pattern  (sometimes  referred 
to  as  fringes)  shown  in  Figure  5.  The  geometry  for  the  array  pattern  formed  by  a  two-element 
interferometer  is  illustrated  in  Figure  6.  The  interference  pattern  can  be  recreated  numerically  if 
the  processing  is  coherent.  This  is  accomplished  if  the  baselines  are  short  and  a  common  local 
frequency  standard  is  used  as  the  reference  for  all  phase  data. 
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TWO-ELEMENT  INTERFEROMETER  ARRAY  PATTERN 
1 0-m  S«(Miratlofl  at  L« 


HGURE  5.  TWO-ELEMENT  INTERFEROMETER  ARRAY  PATTERN 


FIGURE  6.  TWO-ELEMENT  INTERFEROMETER  ANTENNA  ARRAY 
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Given  a  fixed  baseline  defined  by  the  phase  centers  of  two  antennas,  the  array  pattern  is 
fixed  to  and  rotates  with  the  earth.  The  GPS  satellites  will  pass  through  the  fringes  as  they  pass 
overhead.  If  the  signals  from  the  two  antennas  shown  in  Figure  6  were  connected  to  a  normal 
GPS  receiver,  it  could  track  the  satellites  in  the  usual  way  with  the  exception  that  the  antenna 
pattern  has  peaks  and  nulls  rather  than  being  omnidirectional.  If  the  receiver  has  an  accurate  way 
to  record  signal  strength  these  observ’ations  of  signal  strength  for  each  satellite  over  a  period  of 
time  would  be  enough  to  fix  it  in  a  particular  track  through  the  sky.  Since  the  satellite  positions 
are  known,  the  direction  of  the  baseline  determined  in  the  usual  relative  positioning  mode  and 
the  satellite  tracks  in  the  sky  should  be  enough  to  isolate  the  location  of  the  observer.  The 
accuracy  with  which  this  could  be  accomplished  depends  upon  the  null-to-null  spacing  of  the 
array  pattern. 


D^{e)  =  1  +  cos 


'  2  nb  sing\ 

.  ^  I 


(11) 


From  the  equation  for  the  array  pattern  [Equation  (11)],  the  change  Ae  of  elevation  angle 
£  for  one  cycle  of  interference  when  e  <=  W2  is  approximately  k/h.  In  the  ca.se  where  b  is  100  m, 
the  fringe  spacing  is  about  0.002  rad  at  L,.  Figure  5  illustrated  the  fringe  pattern  as  a  function 
of  E  for  a  baseline  of  10  m.  Notice  that  the  spacing  varies  according  to  the  sine  of  e  and  is  coarse 
at  the  horizon  and  fine  at  the  zenith.  The  variation  may  give  some  additional  information  about 
the  elevation  angle  of  the  satellites. 

As  with  the  other  schemes  outlined  above,  the  short  distance  between  antennas  limits  the 
angular  precision  and  consequently  the  accuracy  of  the  derived  ab.solute  position.  In  the  example, 
the  2-mrad  fringe  spacing  that  results  from  a  100-m  baseline  corresponds  to  a  distance  of  40  km 
at  GPS  altitudes.  If  observations  over  a  period  of  time  can  subdivide  the  null  to  null  angle  into 
1000  parts,  the  satellite  position  can  still  be  in  error  by  about  40  m.  This  is  on  the  order  of  the 
error  introduced  by  SA-epsilon  and,  without  further  analysis,  seems  to  offer  little  if  any 
improvement  in  the  absolute  position  one  would  obtain  with  the  conventional  SPS  techniques. 


MININ  TZING  THE  RESIDUAL  OF  FIT 

When  performing  relative  positioning  solutions,  the  satellite  ephemeris  and  an  estimate 
of  the  absolute  position  of  the  reference  antenna  is  required  to  compute  in  Equation  (2).  A 
plane  geometrical  representation  is  presented  in  Figure  7.  For  each  antenna  pair,  the  observation 
is  the  time  difference  of  arrival  (/)  of  the  wavefront  between  the  two  antennas.  This  can  be 
expressed  as  a  function  of  the  baseline  h,  the  range  r,  and  the  elevation  angle  e  of  the  satellite 
as  seen  from  the  reference  antenna.  The  expression  for  /  in  length  units  is  shown  in  Equation 
(12). 
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HGURE  7.  RELATIVE  POSITIONING  GEOMETRY 


-  rbcose  -  r 


(12) 


A  translational  shift  in  the  location  of  the  reference  antenna  or  the  satellite  will  produce 
a  change  in  the  elevation  angle.  The  sensitivity  of  e  to  a  change  in  /  is  given  by  the  partial 
derivative  in  Equation  (13).  For  the  case  of  GPS  positioning  a  point  on  the  earth,  the  range  r>>b 
and  the  partial  can  be  approximated  by  Equation  (14). 


de 


-rbsine 


yj +  b^  -  rb cose 


(13) 
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—  -  -ftsine  (14) 

de  2 


Assuming  the  satellite  position  is  known  through  its  ephemeris,  a  small  horizontal 
displacement  in  the  baseline  As  will  cause  a  change  in  the  elevation  angle  of  approximately  As/r 
(for  e«90  deg).  This  in  turn  leads  to  a  change  in  /  given  by  Equation  (15). 


A/  = 


As 

r 


(15) 


With  e  <=>  90  deg,  r  *»  2x10®  m,  and  b  »  30  m,  the  evaluation  of  Equation  (15)  returns  the 
expression  given  in  Equation  (16). 


A  /  =  7.5x10-^  As 


(16) 


If  the  minimum  detectable  A/  is  1  mm,  then  the  expected  sensitivity  to  an  error  in  absolute 
position  is  about  13.3  km. 


MULTIPLE  BASELINE  EXPERIMENT 

An  experiment  to  investigate  the  potential  of  the  technique  above  was  performed  in 
January  1993  at  NSWCDD.  A  direction-finding  receiver  was  used  with  four  antennas.  The 
particular  receiver  used  had  24  satellite  channels  that  were  allocated  in  four  banks  of  six  L, 
trackers.  The  antennas  were  connected  to  the  receiver  through  30-m  cables.  Thus,  with  the 
receiver  in  the  center  of  a  square,  the  opposite  comers  can  be  up  to  60  m  apart.  Since  the  four 
banks  of  trackers  originate  from  the  same  receiver,  there  is  no  time  or  frequency  deviation 
between  them.  Therefore,  the  standard  double-difference  would  not  be  required;  a  single¬ 
difference  between  satellite  measurements  should  suffice.  A  diagram  of  the  experimental  layout 
is  shown  in  Figure  8. 

Two  consecutive  days  of  data  were  collected  at  the  site.  The  time  interval  between 
observations  was  selected  to  be  20  sec.  The  distribution  of  satellites  on  14  January  1993  was  as 
depicted  in  Figure  9.  The  elevation  angle  below  which  no  data  were  accepted  was  set  to  12  deg. 
The  on-line  receiver  software  selected  and  automatically  tracked  the  same  set  of  satellites  at  each 
antenna.  The  data  were  recorded  internally  and  were  downloaded  once  each  day.  The  standard 
double-difference  software  uses  a  single-reference  satellite  for  each  solution.  When  the  solution 
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is  performed,  the  user  has  the  option  of  choosing  which  satellite  should  be  the  reference,  the 
identity  of  the  other  satellites  included,  and  the  number  of  observation  epochs  to  process. 


Antenna 


Four-Antenna  Aghtech 
3DF  Receiver 


30-m  Antenna 
Cables 


Four-Antenna  GPS 
Positioning  Experiment 
at  NSWCDD 


Building  \2d0 


FIGURE  8.  FOUR-ANTENNA  ARRAY  AT  THE  HERO  SITE  AT  NSWCDD;  14-15  JANUARY  1993 


To  illustrate  the  repeatability  of  the  solutions,  six  different  time  spans  (partially 
overlapping)  with  different  reference  satellites  were  processed  each  day.  The  12  solutions  for 
each  baseline  are  listed  in  the  Tables  1  through  3.  Figures  10  through  12  show  the  scatter  in  the 
solutions  listed  in  the  tables.  One  of  the  trackers  in  the  fourth  bank  of  six  channels  suffered  from 
an  excessive  number  of  unexplained  phase  jumps.  Therefore,  the  satellites  tracked  by  this  channel 
were  omitted  from  the  solutions.  As  a  result,  some  of  the  baselines  were  not  as  well  determined 
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as  they  could  have  been  had  the  data  from  those  satellites  been  processed.  Even  with  the  bad 
satellite  channel,  the  baselines  were  repeatability  determined  at  the  millimeter  level. 


GPS  SATELLITE  VISIBILITY 
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HGURE  9.  SATELLITES  VISIBLE  AT  NSWCDD  ON  14  JANUARY  1993 


Table  1  lists  the  results  for  the  12  solutions  computed  for  baseline  1-2.  The  leftmost 
column  lists  the  day  and  the  second  column  shows  the  reference  satellite  used  for  that  solution. 
The  next  three  columns  are  the  x,  y,  z  components  of  the  baseline  vector  followed  by  its  length. 
The  rightmost  column  lists  the  residual  for  the  double-difference  fit.  The  average  residual  was 
found  to  be  3.010  mm  with  a  standard  deviation  of  0.220  mm.  The  residual  was  a  little  larger 
than  the  1.9-mm  standard  deviation  on  the  baseline  length.  Tables  2  and  3  show  similar  results 
for  baselines  1-3  and  1-4,  respectively.  The  average  residual  for  baseline  1-3  was  3.088  mm  and 
for  baseline  1-4  it  was  3.083  mm. 
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The  root-mean-square  (rms)  of  the  residuals,  however,  were  consistently  between  2  and 
4  mm.  The  residual  plots  show  that  this  rather  high  rms  is  due  to  systematic  variations  that  are 
not  attributable  to  system  noise.  These  signals  may  be  due  to  multipath  effects,  some 
environmental  condition,  or  a  propagation  effect  that  repeats  from  day  to  day.  Examples  of  the 
summary  residual  plots  for  each  baseline  and  each  satellite  pair  on  day  015,  with  PRN  19  as  the 
reference  satellite,  are  shown  in  Figures  12  through  15.  A  comparison  from  one  day  to  the  next 
for  one  satellite  pair  (PRN  11-19)  is  illu.strated  in  Figures  16  and  17.  Notice  how  well  this 
residual  is  repeated  on  the  second  day;  this  signature  is  typical  of  multipath.  However,  it  does 
not  explain  why  the  other  satellite’s  residuals  in  the  summary  plots  are  so  similar  to  the  PRN  1 1- 
19  residual.  Though  all  the  residuals  in  Figure  14  have  PRN  19  in  common,  it  is  surprising  that 
the  similarity  is  so  strong.  Similar  comments  can  be  made  about  the  baseline-to-baseline 
correlations. 

For  this  demonstration,  repeated  relative  positioning  solutions  were  performed  keeping  all 
variables  the  same  except  for  the  given  site  position.  This  mimics  having  the  site  position 
unknown  and  letting  the  rms  of  the  relative  positioning  solution  be  the  indicator  of  the  true 
position.  The  lower  the  residual  of  the  fit,  the  smaller  the  model  error  due  to  position  error. 
Given  many  starting  positions  (each  one  offset  a  different  amount  from  the  correct  value),  the 
rms  of  the  residual  should  indicate  a  minimum  value  near  the  true  location.  This  is  verified  by 
the  results  for  longitude  shown  in  Figure  18  and  for  latitude  shown  in  Figure  19.  In  each  case, 
the  minimum  is  too  shallow  to  accurately  determine  the  true  location  to  better  than  about  1  km. 
Accuracies  on  the  order  of  10  m  would  require  lengthening  the  baselines  by  a  factor  of  100.  The 
30  m  baselines  would  need  to  become  3  km  or  greater.  Short  baselines  have  the  advantage  that 
real-time  solutions  could  be  obtained  on  the  spot.  When  receivers  are  separated  by  a  few 
kilometers,  a  communication  link  would  need  to  be  employed.  Near  real-time  precise  relative 
positioning  software  has  become  available  for  baselines  up  to  50  km. 
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TABLE  1.  SOLUTIONS  FOR  BASELINE  1  -  2 


Day 

Reference 

Satellite 

X(m) 

Y(m) 

Z(m) 

L{m) 

Residual 

(mm) 

014 

13 

20.8746 

30.1871 

31.5376 

48.3904 

2.656 

014 

19 

20.8740 

30.1881 

31.5383 

48.3912 

2.866 

014 

23 

20.8735 

30.1889 

31.5375 

48.3909 

2.941 

014 

24 

20.8745 

30.1876 

31.5375 

48.3906 

2.878 

014 

25 

20.8746 

30.1874 

31.5380 

48.3909 

3.071 

014 

32 

20.8734 

30.1885 

31.5365 

48.3900 

3.032 

015 

13 

20.8741 

30.1855 

31.5385 

48.3898 

3.232 

015 

19 

20.8734 

30.1875 

31.5391 

48.3911 

2.890 

015 

23 

20.3740 

30.1885 

31.5385 

48.3915 

2.730 

015 

24 

20.8740 

30.1856 

31.5394 

48.3904 

3.218 

015 

25 

20.8742 

30.1867 

31.5395 

48.3912 

3.278 

015 

32 

20.8736 

30.1869 

31.5392 

48.3909 

3.326 

Average 

20.8740 

30.1874 

31.5383 

48.3907 

3.010 

Standard 

Deviation 

0.0004 

0.0011 

0.0009 

0.0005 

0.220 

FIGURE  10.  RESULTS  OF  TABLE  1  PLOTTED 
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Average 


Standard 

Deviation 


Reference 

Satellite 


TABLE  2.  SOLUTIONS  FOR  BaSELINE  1  -  3 


X{m) 


48.8642 


48.8640 


48.8636 


48.8638 


48.8642 


48.8635 


48.8637 


48.8632 


48.8629 


48.8630 


48.8637 


48.8629 


48.8636 


0.0005 


Y(m) 

Z(m) 

L(m) 

25.8395 

18.1896 

58.1915 

25.8406 

18.1898 

58.1919 

25.8413 

18.1903 

58.1920 

25.8400 

181895 

58.1914 

25.8396 

18.1895 

58.1915 

25.8411 

18.1895 

58.1916 

25.8374 

18.1913 

58.1907 

25.8403 

18.1909 

58.1914 

25.8412 

18.1906 

58.1915 

25.8378 

18.1916 

58.1904 

25.8395 

18.1917 

58.1918 

25.8396 

18.1917 

58.1911 

25.8398 

18.1905 

58.1914 

0.0012 

0.0009 

0.0005 

Residual 

(mm) 


nOURE  11.  RESULTS  OF  TABLE  2  PLOTTED 
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Day 

Fteference 

Satellite 

X(m) 

1 

Y(m) 

Z(m) 

L(m) 

Residual  (mm) 

014 

13 

30.8107 

4.3793 

-14.1850 

34.2007 

2.607 

014 

19 

30.8110 

4.3746 

-14.1859 

34.2008 

4.145 

014 

23 

30.8096 

4.3763 

-14.1865 

34.2000 

2.817 

014 

24 

30.8107 

4.3798 

-14.1845 

34 '’007 

2.734 

014 

25 

30.8090 

4.3774 

-14.1821 

34.1978 

3.770 

014 

32 

30.8093 

4.3774 

-14.1869 

34.2000 

3.184 

015 

13 

30.8098 

4.3803 

-14.1835 

34.1995 

2.870 

015 

19 

30.8088 

4.3779 

-14.1850 

34.1989 

2.818 

015 

23 

30.8099 

4.3773 

-14.1852 

34.1998 

2.760 

015 

24 

30.8098 

4.3804 

-14.1831 

34.1993 

2.904 

015 

25 

30.8006 

4.3803 

-14.1835 

34.1993 

3.329 

015 

32 

30.8094 

4.3789 

-14.1841 

34.1992 

3.057 

Average 


Standard 

Deviation 


30.8098 

-4.3783 

-14.1846 

34.1997 

3.083 

0.0007 

0.0018 

0.0014 

0.0009 

0.461 

RGURE  12.  RESULTS  OF  TABLE  3  PLOTTED 
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HGURE  13.  SUMMARY  PLOT  FOR  BASELINE  1-2  ON  DAY  015  WITH 
PRN  19  AS  REFERENCE 


nOURE  14.  SUMMARY  PLOT  FOR  BASELINE  1-3  ON  DAY  015  WITH 
PRN  19  AS  REFERENCE 
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nOURE  15.  SUMMARY  PLOT  FOR  BASELINE  1-4  ON  DAY  015  WITH 
PRN  19  AS  REFERENCE 


nOURE  16.  SATELLITE  PAIR  11-19  DOUBLE-DIFFERENCE  RESIDUAL 
FROM  BASELINE  1-3  ON  DAY  014 
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HGURE  17.  SATELLITE  PAIR  11-19  DOUBLE-DIFFERENCE  RESIDUAL 
FROM  BASELINE  1-3  ON  DAY  015 


FIGURE  18.  DOUBLE-DIFFERENCE  RMS  RESIDUAL  FOR 
LONGITUDE  OFFSETS 
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OPPS<AUTO)  HESIDUAl.  ERPOP  VS  ASSUMED  POSITION  ERROR 


HGURE  19.  DOUBLE  DIFFERENCE  RMS  RESIDUAL  FOR  LATITUDE 

OFFSETS 


GPS  AND  LEVELING 

The  objective  of  this  procedure  is  to  obtain  a  vector  normal  to  the  ellipsoid  at  the 
reference  site.  This  vector  defines  the  latitude  and  longitude  at  the  site.  The  procedure  to  estimate 
the  vector  components  is  given  below  and  is  demonstrated  using  data  collected  during  a  test  at 
NSWCDD  in  1987. 

The  procedure  involves  collecting  both  GPS  and  first-order  leveling  data  along  two 
baselines  radiating  from  the  reference  site  and,  preferably,  separated  by  an  angle  of  about  90  deg. 
Using  an  assumed  initial  value  of  position,  the  baseline  vector  (defining  Ax,  Ay,  and  Az)  is 
determined  using  standard  GPS  relative  positioning  techniques.  These  vector  values  have  been 
found  to  be  somewhat  insensitive  to  the  initially  assumed  position  for  the  baseline  lengths 
discussed  here.  Use  of  the  C/A-determined  initial  position  is  adequate.  The  baseline  length, 
heights  (h,  and  /i^)  above  the  ellipsoid,  and  azimuth  are  then  obtained.  In  addition,  first-order 
leveling  difference  in  height  is  obtained  (AH)  and  values  of  deflection  of  the  vertical  (^  and  i^) 
are  obtained  using  a  gravitional  model.  The  deflection  in  the  baseline  direction  is^ 
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e  =  (cosff)  (  +  (sinff )  rj 


(17) 


where 

a  =  geodetic  azimuth 

I  =  north-south  component  of  vertical  deflection 
r]  =  east-west  component  of  vertical  deflection 

The  height  correction  to  ^or  the  deflection  e  is  given  by  Equation  (18),  where  is  the 
baseline  length. 

A/i^  =  -eAs  (1^) 


The  leveling  AH  and  vertical  deflections  Ah^  corrections  are  made  to  the  height  If  the 
geoid  and  the  ellipsoid  model  did  not  curve,  these  would  be  the  only  corrections  needed  to  obtain 
the  tangent  vector  to  the  ellipsoid.  However,  a  curvature  correction  is  required  as  illustrated  in 
Figure  20.  An  approximate  geometric  correction  is  given  by  the  following  expressions  where  7?^ 
is  the  radius  of  the  earth. 


^  Re 


Ah 


cosd 


(19) 


'  K  -  -  A*.  ■* 

The  position  of  the  tangent  vector  can  now  be  described  by  Equation  (20)  and  the  original 
latitude  and  longitude  by  <j»2,  These  values  are  converted  to  earth-fixed  orthogonal  coordinates 
of  X2,  y^,  and  Z2  and  differenced  with  the  reference  site  values  to  obtain  the  ellipsoid  tangent 
vector  defined  as  v^,  v^;,  and  v^.  Next,  a  second  tangent  vector  is  obtained  by  the  above 
procedure  for  a  baseline  from  the  reference  site  and,  preferably,  roughly  perpendicular  to  the  first 
baseline.  The  second  vector  is  def.ned  as  and  v^.  The  cross  product  is  given  by  the  vector 

w  whose  components  are  expander  in  Equations  (21)  and  (22).  The  corresponding  unit  vector 
components  are  found  using  Equation  (22). 
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The  last  step  of  the  procedure  is  to  obtain  the  estimated  latitude  and  longitude  (^,  and 
from  the  normal  unit  vector  relationship  given  in  Equation  (7).  Explicitly  this  is  given  in 
Equation  (23). 


=  arcsin(u3) 

f  i  “2 

X^  =  arctan  — 


(23) 


The  positioning  procedure  has  been  demonstrated  and  validated  using  test  data  taken  in 
1987  at  NSWCDD.®  The  GPS  baselines  and  the  leveling  loop  are  presented  in  Figure  21.^  For 
the  current  demonstration,  only  two  baselines  were  analyzed;  these  are  from  site  MERE  to  BOM2 
and  MERE  to  CHUR.  At  the  reference  site  (MERE),  astrogeodetic  vertical  deflection 
measurements  were  also  obtained  by  the  U  S  Naval  Observatory  (USNO).  Two  and  three  days 
of  GPS  data  were  collected  at  the  above  baselines,  respectively,  using  TI4100  receivers.  These 
data  were  processed  using  the  National  Geodetic  Survey  (NGS)  PHASER  software  and 
DMA/NSWCDD  post-fit  precise  ephemerides.  The  vector  values,  the  baseline  length,  and  the 
azimuth  from  the  GPS  solutions  are  given  in  Table  4;  also  given  are  the  first-order  leveling 
height  differences. 


TABLE  4.  RESULTS  FROM  THE  TEST  AT  NSWCDD 


Vector  Values 

BASELINE 

MBRE  to  BOM2 

MBRE  to  CHUR 

Ax  (m) 

464.891 

-4668.231 

Ay  (m) 

2318.125 

-368.816 

Az  (m) 

2720.982 

913.228 

As  (m) 

3640.660 

4770.995 

Aa  (m) 

15.667219 

283.866708 

AH  (m) 

-1.8364 

27.7834 

The  demonstration  analyses  consists  of  two  cases.  The  normal  to  the  ellipsoid  vector  is 
found  with  and  without  a  qualified  position  offset.  The  position  offset  must  be  qualified,  since 
it  was  not  used  for  both  the  solutions  of  the  baseline  vector  Ax,  Ay,  and  Az  or  for  the  modeled 
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vertical  deflection  values;  these  are  both  insensitive  to  assumed  position.  The  modeled  vertical 
deflection  values  =  -2.89  arcsec  and  t]  =  4.79  arcsec)  were  obtained  from  the  Rapp  360x360 
set  of  potential  coefficients  model  determined  from  measured  gravity  anomaly  data.* 

The  demonstration  results  are  given  in  Table  5.  The  resulting  position  estimate  error  is 
17.8  m  in  latitude  and  21.7  m  in  longitude.  These  errors  appear  to  be  insensitive  to  the  initial 
position  offset  of  100  m  in  each  of  the  x,  y,  and  z  component  directions.  Computations  were 
performed  for  several  additional  cases  as  a  means  to  quantify  the  errors  in  the  above  estimates. 
These  analyses  assume  the  true  location  as  an  initial  starting  values  as  in  case  2  of  Table  4.  The 
first  correction  is  to  use  measured  values  for  vertical  deflections.  These  values  (|  =  -2.52  arcsec 
and  Ti  =  5.30  arcsec)  were  obtained  by  the  USNO’  using  the  Danjon  astrolabe  to  an  accuracy 
of  about  0.1  arcsec.  This  improves  the  case  2  results  to  estimated  position  errors  of  6.3  and  5.9  m 
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in  latitude  and  longitude,  respectively.  Consequently,  the  largest  error  sources  in  the  computation 
are  the  modeled  vertical  deflection  errors.  Next,  case  2  is  improved  further  by  replacing  the 
combined  leveling  and  vertical  deflection  height  difference  values  by  the  difference  in  height 
from  GPS  values  (Note:  this  could  not  be  done  at  the  offset  point,  since  is  very  dependent 
on  the  assumed  position).  This  reduces  the  estimated  position  error  to  3.7  m  and  3.6  m  in  latitude 
and  longitude,  respectively.  This  indicates  that  leveling  is  the  second  largest  error  source  in  the 
procedure  and  that  GPS  contributed  relatively  small  errors. 


TABLE  5.  POSITION  ERRORS  OF  PROCEDURE  DEMONSTRATIONS  USING  THE  TEST  DATA 


Case 

Assumed 
Position 
Offset*  (m) 

Vertical 

Deflection 

Latitude  Error 
(m) 

Longitude 
Error  (m) 

1 

100 

Modeled 

17.798 

21.660 

2 

0 

Modeled 

17.798 

21.660 

*  See  discussion  in  text. 


SUMMARY 

Several  techniques  that  have  the  potential  to  determine  absolute  positions  in  the  presence 
of  SA  have  been  considered.  All  use  some  form  of  differencing  of  the  GPS  phase  data  to  remove 
the  dither  component  from  SA.  Tho.se  that  require  an  auxiliary  measurement  are  less  sensitive 
to  the  epsilon  component  also.  The  two  primary  limiting  factors  that  prevent  these  techniques 
from  being  practical  are  the  precision  with  which  the  time  difference  of  arrival  (TDOA)  of  the 
wavefront  at  the  ends  of  the  baseline  can  be  measured.  The  current  limit  is  about  0.005  cycle  or 
1  mm  at  L,.  The  other  limiting  factor  is  the  distance  between  the  antennas.  Since  the  TDOA 
measurement  is  essentially  independent  of  the  baseline  length,  a  larger  antenna  separation  gives 
a  correspondingly  more  precise  angular  measurement.  Baselines  on  the  kilometer  scale  or  larger 
are  required  for  accurate  estimates  of  the  absolute  position. 
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